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Abstract 

Fatigue crack advance induced by the application of cyclic quasistatic loads is 
investigated both numerically and analytically using a lattice spring model. 
The system has a quasi-one- dimensional geometry, and consists in two sym- 
metrical chains that are pulled apart, thus breaking springs which connect 
them, and producing the advance of a crack. Quasistatic crack advance oc- 
curs as a consequence of the plasticity included in the springs which form 
the chains, and that implies a history dependent stress-strain curve for each 
spring. The continuous limit of the model allows a detailed analytical treat- 
ment that gives physical insight of the propagation mechanism. This simple 
model captures key features that cause well known phenomenology in fatigue 
crack propagation, in particular a Paris-like law of crack advance under cyclic 
loading, and the overload retardation effect. 
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1. Introduction 

Crack propagation usual ly occur s in o ne of two qualitatively different 



forms. In dynamic fracture ( iFreundl . ll99Cf ). propagation is very rapid, typ 



ically with velocities that are an appreciable fraction of the sound velocity 
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in the material. Dynamic propagation occurs when some threshold load 
is exceeded. When this occurs, it becomes energetically favorable for the 
crack to advance. Once propagating, dynamical effects limit the maximum 
attainable crack velocity to some fraction of the sound velocity. Veloci- 
ties in the range of 30 - to-60 % of the sound velocity are typically observed 
(IFineberg and Marderl . Il999l ). 

Even if the threshold load for dynamical propagation is not reached, there 
is still the possibility of what is called sub-critical crack growth. This kind 
of slow propagation can occur through different mechanisms. When the 
external load is kept constant, two different mechanism may lead to crack 
growth: creep and stress corrosion. Creep occurs at elevated temperatures, 
by degradation of the material ahead of the crack tip and corrosion occurs 
by a material-environment interaction close to the crack tip. Both of them 
are strongly time dependent. 

Another class of sub-critical crack growth, is fatigue. More precisely, 
cyclic fatigue, which is the main focus of this paper. It occurs when the 
external load has a cyclic dependence in time. We limit the analysis to the 
cases where the external load changes very slowly in time, and so dynamical 
effects can be ignored altogether. We also assume that a seed crack already 
exists in the system at the beginning of the process, as we do not intend 
to analize the problem of crack nucleation. A necessary condition for an 
applied cyclic stress to produce crack advance is to reach conditions in which 
plastic yielding in the material occurs, at least in some neighborhood of 
the crack tip. Plastic effects make stress-strain curves of the material be 
history dependent. So, after a cyclic variation of the external load, the state 
of the system may not be identical to that at the beginning of the cycle. 
Particularly, a finite crack advance in each cycle can occur. It is important 
to emphasize that in this case the process does not require the existence of 
any activation step, allowing a fully deterministic study of the problem once 
the constitutive elasto-plastic behavior of the material is known. 

In general terms, for a mechanical component, a stage of cyclic fatigue 
propagation can eventually lead to a regime of dynamic propagation, at which 
abrupt failure occurs. Cyclic fatigue propagation may be the main concern 
in applications in which mechanical components are exposed to temporally 
variable loads, or to repetitive thermal cycling. 

We focus on cyclic fatigue propagation in the present paper. Although the 
amount of experimenta l data and the phenomenolo gical treatments of this 
problem are abundant ( ISureshl . Il998l ; iBolotinl . Il999l ) , the basic studies that 
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are available on cyclic fatigue crack propagation are scarce. The difficulty 
of the problem is the typical one of fracture mechanics (i.e., the spanning 
of many orders of magnitude between the process zone at the crack tip to 
the macroscopic region described by the continuum dynamics) plus the ne- 
cessity of a det ailed descript i on of plasticity in the material. For instance, in 



a recent work (IFarkas et al.l . 120051 ) , crack advance under cyclic loading has 



been obtained in atomistic simulations of a nano-crystalline material. The 
plasticity in this case is seen to be related to dislocation emission from the 
crack tip upon loading. Results of this kinds of studies are at present lim- 
ited by the availability of computing resources. An alternative approach is 
to consider the system as des cribed by the equa tions of continuum media 



that include plastic response ( iMoes et al.l . Il999h . However, this approach 



necessarily breaks down sufficiently close to the crack tip, and has to be 
complemented either by full atomistic simulations in a small neighborhood 
of the tip, or by phenomenological prescriptions about the crack tip advance 
behavior. 

From a conceptual point of view, it is always desirable to have simple 
models that reproduce the available phenomenology with a minimum of in- 
gredients. Having in mind the cyclic loading fatigue crack advance prob- 
lem, the possibility that we have been studying corresponds to spring lat- 
tice systems with plasticity in the springs. In fact, lattice spring models 
have been a very important benchmark where many predictions of fracture 



analytic treatments were obtained ( 


Slepvan. 


1981: 


Marder and Gross. 


1995: 


Kessler and Levine. 


1999b. 


2001: 


Guozden and Jaela. 


2006 


). These studies 



have mainly focused on propagation in lattices with linear, non-linear elastic, 
or visco-elastic springs. In all these cases there is a unique relation between 
a stationary applied strain and the stress produced. To our knowledge there 
have not been previous attempts to study cyclic fatigue crack growth using 
lattice spring models with plasticity. 

In the present paper we study the simplest case of a lattice of elasto- 
plastic springs that may cause fatigue crack growth. This is a quasi-one- 
dimensional model consisting of two identical, parallel elasto-plastic chains, 
joined by breakable springs, that are laterally pulled apart by the external 
load, in such a way that a crack can propagate between chains. Formally, 
this represents a case of "Mode III" propagation, as we only consider per- 
pendicular displacements to the chain direction. In absence of plasticity, this 
kind of "one-dimensional" geometry has been studied in detail in the context 
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of dynamical fracture, and has proved to be useful as a simple benchmark 
for the more complex beha vior that is ob s erved in a more thorough t wo di- 
mensional implementation ( iLangerl . Il992t iBouchbinder and Lol. 



have also used this model in a previous work f lGuozden and Jaglal . 120061 ) and 



2008). We 



obtained results in the dynamical propagation case. We now concentrate in 
this quasi one-dimensional model because in addition to the numerical imple- 
mentation, it allows a detailed, mostly analytical description of the fatigue 
propagation mechanisms involved, and also because in this case, the germ of 
experimental features of fatigue crack growth is already observed. We leave 
for a forthcoming publication the study of a more realistic two dimensional 
geometry, which can be made only through numerical simulation. 

In the next Section we describe the quasi-one-dimensional model and the 
numerical technique in detail. In Section III we present the main results that 
we have obtained. They include the very observation of the possibility of 
fatigue crack advance under cyclic conditions following a Paris-like law, and 
the observation of the overload retardation effect, in which a cycle using a 
deformation larger than the average induces an eventual retardation in the 
advance of the crack. Finally in Section IV we discuss the results of our 
approach in the context of other studies and present the final conclusions. 



2. Details of the model and the numerical technique 

The system we model is reminiscent of a stripe geometry, under mode III 
conditions imposed by rigid displacements on the lateral sides, with a crack 
advancing in the middle of the stripe. In this quasi one- dimensional idealiza- 
tion (see Fig. [Q we consider only two symmetric chains around the middle 
line of the system: u(x) and —u(x). They are respectively connected through 
linear springs to the lateral strain gauge located at ±35/2. In addition, the 
two middle chains are connected by breakable springs. As soon as the length 
of any of these springs exceeds some threshold value u^k, the spring breaks 
irreversibly. We emphasize that these are the only springs allowed to break 
in the model. The intact inter-chain springs are shadowed in Fig. [TJ 

The equations of the model are those of mechanical equilibrium for a 
given profile of the chain u(x). To simplify, we first consider the case in 
which the continuous profile u(x) is replaced by a discrete set of values Uj, 
where j corresponds to a horizontal coordinate j A, A being the discretization 
parameter. The equilibrium equation for Uj reads 
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Q<5 - Uj \ A - 26v,jA + F(u j+1 - uj) + F{uj- X - uj) = 0. (1) 

The different terms represent: 

1) the coupling of Uj to the upper border located at 35/2. 

2) The coupling of Uj to the mirror chain at —u(j). In this term the 
Heaviside function 6 is defined as one, unless 2uj was larger than at some 
previous time, in which case it is taken as zero. Note that in the initial 
condition we take 6 = for all j lower than some value, to simulate a pre- 
existent crack. The factor A included in terms 1) and 2) corresponds to use 
a spring constant that is unitary for a unit length of the continuous system. 

3,4) the last two terms are the intra-chain forces with the particles at the 
left and at the right of the site j. Plasticity is included precisely in these 
terms, and not in the interchain springs. This choice looks strange at first, 
since springs that eventually break are those which experience the largest 
deformation, and plastic deformation is expected to be maximum for them. 
In spite of this, we have found that the main characteristics of cyclic fatigue 
can be discussed disregarding plasticity of inter chain springs. In other words, 
its inclusion does not modify the results in a qualitative way. Since this is 
the key simplification that allows a detailed analytical treatment, we think 
that this advantage justifies neglecting plasticity in vertical springs. 
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Figure 1: A sketch of our idealized system. The crack tip position x = Xq separates regions 
I and II. Intact inter-chain springs are shadowed. Note that h is given by h= **f^-. The 
part of the chains that has accumulated some plastic deformation is shown by a thicker 
line (note that the region with plasticity penetrates slightly into region I) . 

We now concentrate in the intra-chain forces. For its description, we split 
the force F in two parts, an elastic part f e and an inelastic part f p where the 
effects of plasticity will be included. The total force will be F = (1 — P)f e + 
Pf p where the parameter P (0 < P < 1) controls the extent of plasticity 
of the spring. The elastic part is simply calculated as f e = (uj±i — Uj)/A. 
The spring constant is rescaled with the discretization in such a way that it 
would be one if A = 1. To specify the inelastic part, we introduce a "rest 
length" L°jj + \ for the spring connecting sites j and j + 1. We calculate the 
inelastic force as f p = (uj — Uj + i — L° J+1 )/A. The total force in the spring 
is thus given by 

F=[(l- P)(u j+1 - Uj ) + P(u 3 - u j+1 - L° j>j+1 )] /A (2) 

The force f p models an ideal plastic behavior. The rest length L is 
initially set to zero to model a virgin sample. Upon small variations of 
Uj + i — Uj, L is changed, if necessary, to avoid f p to go outside some pre- 
established range ±u n iA, namely, if \f p \ becomes larger than n n /A, L Q is 
adjusted to get \f p \ = u n iA. We will refer to cases where values of L different 
from zero appear, as cases in which a "plastic deformation" is present in the 
system. 
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Through this mechanism, we can see that history dependent forces appear 
in the system. In Fig. ([2]) we show the evolution of f p and the total force F 
as the length of the spring / = Uj+i — Uj is changed in a prescribed manner. 
Note that although the evolution is fully deterministic, the force is not a 
single valued function of I. For instance, the spring has the same length I at 
points 1 and 4 but the force it exerts is different. 

The numerical procedure we use to find stationary solutions to this prob- 
lem consists in starting with a pre-existent crack and no plastic deformation 
(Lq = 0) in the system, and increase the value of 5 in small steps. At each 
step the solution to the equilibrium equations is obtained through a relax- 
ation protocol, and after that the value of 5 is increased again. Once some 
maximum value of 5 is reached, the same procedure is repeated for decreasing 
5. As it was already emphasized, the history dependent stress-strain relation 
of the springs makes it possible to observe a systematic crack tip advance 
upon cyclic loading conditions. 

The previous description is appropriate for numerical implementation, 
but for an analytical treatment its is better to return to the continuous 
description. This is achieved by letting the discretization parameter A go to 
0. The continuous equation for our model is obtained as 

d 2 u f 35 \ dl 



2 ^ + P~£ = 0, (3) 



dx 2 \ 2 J dx 

where we have defined Iq = Lq/A. In this description, the system is con- 
tinuous along the chains, although it remains essentially discrete in the per- 
pendicular direction. Th is kind of continuous lim it along a single spatial 



dimension is well defined (IGuozden and Jaglal . 120061 ) , contrary to the case of 



a full continuous limit in which the transition from discrete to continuous 
is much more subtle. Given a distribution of plasticity lo(x) the previous 
equation can be solved. However, upon a variation of 5, this equation does 
not stand alone, as it has to be complemented with the prescription that Iq 
adjust itself to satisfy \du/dx\ — Iq(x) < u n[ . 



3. Results 

3.1. Crack advance upon monotonous load increase 

In the absence of plasticity (P = 0), the system is perfectly elastic, and 
an energy balance analysis can be applied: if the strain 5 is lower than some 
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a) 




Figure 2: Qualitative behavior of the plastic springs. In a) the inelastic force f p of a spring 
when / first increases and then decreases is shown. In b) the corresponding evolution of 
length (/), rest length (Iq = Lq/A), and f p , along this path is indicated. Note that when 
the f p reaches its maximum allowed value, Iq starts to change. In c), the total force, 
including the clastic part, is indicated (we use P = 0.5). 
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critical value 5q, there is only a static solution in which the crack tip position 
is stable (assuming crack healing does not occur), as the external load is not 
able to provide the necessary energy for crack advance. Instead, for 5 > 5g, 
there is in principle enough available energy for the system to be in a run 
away state, in which the crack moves forward at some finite velocity. This 
corresponds to a continuous description. The discrete ness of the syst e m in- 



troduces the additi o nal ingredient of lat t ice trapping ([Thomson et all. 11971 



Paskin et al.. 1981; Kessler and Levine. 1999a 



Bernstein and Hessl . 20031 ). 



and the minimum necessary value of 5 to have a run away solution becomes 
somewhat larger than the value 5g determined by purely energetic arguments. 
When plasticity is incorporated into the springs (P > 0), a new regime oc- 
curs, intervening between the static and dynamic regimes. In fact, at a first 
critical value of 5 that we call S G , the crack starts to elongate. However this 
elongation is not unstable now. Instead, there is a well defined value of the 
crack tip advance X as a function of 5. The crack only runs away when a 
second critical value Sq is exceeded. The numerical results that show this 
behavior are presented in Fig. [3J For any value of S between 5 G and 5q the 
crack tip position is stable. 



5 | — 

4 - 
3 

2 - 

1 - 

— 
2.3 



S'n 




2.35 



2.4 

<5/"nI 



2.45 



2.5 



Figure 3: Crack tip position X as a function of an increasing value of 5, starting from a 
virgin configuration (namely, no plastic deformation in the initial state). Parameters are 
u bk /u n i = 4, P = 0.5, and A = 1/20. 

The actual configurations of the system at increasing values of S are pre- 
sented in Fig. HI There we plot the chain profile u(x), and also the cor- 
responding values of the rest length Iq(x) of the springs forming the chain. 
Upon crack tip advance, the portion of the chain that is left behind remains 
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with a finite plastic deformation (i.e. Iq 7^ behind the crack tip). This is of 
the utmost importance as it implies the necessity of some energy expenditure. 
In fact, Fig. [5] shows the typical variation of the plastic force on each spring 
after the crack tip passed. The area of this curve is the energy dissipated in 
the system due to plastic deformation. This is the reason why the crack does 
not immediately destabilize when 5 exceeds the first threshold 5q. 
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Figure 4: Upper curves: chain profiles as 5 increases. Lower curves: corresponding plastic 
deformation Iq of the springs. The thin vertical lines indicate the crack tip position. Note 
that for large values of crack tip advance, the value of Iq behind the crack tip reaches an 
asymptotic value (parameters as in Fig. [3]). 
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Figure 5: Inelastic force on the springs, as the crack passes by. The shadowed area 
corresponds to the energy used in the plastic deformation, that amounts to P(lau n i +Z§/2). 

We proceed now to a detailed analytical study that gives insight into the 
process of quasistatic crack advance upon load increase. In order to obtain 
analytical expressions, we have to consider the continuous limit, namely, the 
case A 0. First of all, we calculate S G and Sq in this limit. 

The calculation of S G can be made as if springs were non-linear elastic, 
since plasticity is not involved. In fact, upon increasing S, as long as no 
spring is broken (as it happens below S G ), all spring lengths also increase. 
This means that springs explore only the ascending part in Fig. [2]2, which 
corresponds to a non-linear elastic behavior. In these conditions, the thresh- 
old can be determined by an energy balance argument, equating the available 
energy ahead of the crack with the energy that is necessary to break a spring. 
This gives the expression for 8 G in our geometry 

3 m = < (4) 

from which 

Sh = ^ (5) 

Note that this value does not depend on the parameter P. 

As soon as the crack tip moves on, some springs reduce their length, ex- 
ploring the descending part of the curve in Fig. [2^, and plasticity becomes 
involved. The calculation of Sq can still be made through an energy bal- 
ance argument, if we take into account the additional energy spent in the 
process. For each individual spring within the chain, its evolution during 
the whole process implies an initial stretching, and then a reversion to a 
state of zero length, as indicated in Fig. |5j The energy Eq spent in this 
process is proportional to the shadowed area in the figure, and its value is 
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Eq = P(l u n i + Iq/2). Knowing the value of l , 5q can be calculated from 
the energy balance equation 



my 



II 



Ml 

2 



2E (Z ) 



(6) 



(the factor 2 in the last term comes from the fact that energy is spent in the 
plasticity of springs in the two mirror chains). 

The value of Iq can be calculated through the following argument: as 
5 — > 8q every spring behind the crack tip gets a plastic deformation Iq 0. 
This means that this part of the chain is effectively linear elastic. Its profile 
is given by 

35 ( 35 u bk \ x 



ulx) 



where the slope of the chain at the crack tip (located at xq = 0) is du/dx 
(ubk — 35) /2. Then, the plastic deformation at the crack tip is, in norm, 



35 



Ubk 



Unl- 



it) 



Inserting this result into Eq. ^ for 5 — > 5q we get 
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(9) 



It can thus be seen that 5q coincides with 5q for P = 0, and becomes 
progressively larger when P increases, allowing a finite range 5q — 5q of 
stable crack elongation. For the particular value P = 0.5 used in previous 
figures the value of 5q is 



U 



nl 



(10) 



Now we derive an equation from which the full curve X(5) can be ob- 
tained. We will refer to the definitions in Fig. CD In particular we consider 
region I as that ahead of the crack tip (x > xq), and region II as that behind 
the crack tip (x < xq). 



1 This is true if ^ < 2 + \/12 — 9P. We restrict to this case throughout this work. 
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Region I: Here the horizontal springs deform monotonically, so again 
plasticity does not take part. We only have to consider that springs behave 
as piecewise linear, with two different spring constants, for elongation below 
and above u n \. This offers the possibility to obtain a solution in this region, 
and to evaluate the slope of the chain at the crack tip. H After some algebra 
we obtain 



dui 
dx 




R(S), (11) 



where we have renamed the expression to R(8) for future use. 

In region II we proceed as following: As we indicated previously in the 
introduction, the existence of a non-zero Iq(x) function can be incorporated 
into the equilibrium equation of the system as and additional force Pdl /dx. 
This means that the equation describing the deformation u(x) in region II is 

llt . f 35 , A T-,dlo 



u »( x )+^--u(x))+P^, (12) 
with boundary conditions 

U(-OO) = y, U(X ) = (13) 

Solving formally this equation we find the derivative of the chain at the 
crack tip in terms of Iq(x). We obtain 



du 



n 



dx 



X=XQ J ~°° 



p I X ° %-^o) dy _ ^±^k (14) 

dy 2 



It is interesting to emphasize the meaning of this equation, stating that any 
non-homogeneous Iq(x) behind the crack tip contributes to the slope of the 
chain at the crack tip, and its effect decays exponentially with distance. 
Equating the results for regions I (Eq. [TT]) and II (Eq. [T4"j) we obtain 

P r ^e-^dy - = R(5). (15) 

j-oo dy 2 



2 Here we assume that the chain actually explores the nonlinear regime. This happens 
when |2i^V3| > u n i 
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To solve this integral equation, we multiply each of the terms by e x ° and 
then take a derivative upon a variation in the crack tip position xq. In doing 
this, we also use the fact that upon an advance of the crack tip, the values 
of Iq are not modified behind the crack tip. [f| The result is 



dl 

-r- 1 

axo 



3 d5 
2 dxo 



Z8-u hk 
2 



dR(5) 
dxo 



+ R(S) 



(16) 



The value of Iq at the crack tip in the first term Iq(xo) can be written as 
lo(x ) = du/dx 
we obtain 



X=XQ 

dR{5) 



P 



— u n i = R(5) — u n i and after replacing in equation (JIB]) 
3 d5 35 — Ubk dR(S) 



+ R(5) 



(17) 



dxo 2 dxo 2 dxo 

Changing from the strain variable 5 to h, as defined in Fig. flTJ, we arrive to 



t dR{h) 

dXn 



dh 

3— + 3h 

dx 



Ubk 



dR(h) 
dx 



R(h), 



where R{h) is now 



R(h) 



^3h 2 - Pu 



In addition, using Eq. ( ITTjl . the derivative of R with respect to x is 



dR(h) 



h— 

dx 



dx l-PR(h) 
and combining with Eq. ffTS]) we obtain 



(18) 



(19) 



(20) 



dh 
dx 



R{h) 
h 



3-3 



-h + 



Ubk 



R(h) 



3-3 



(21) 



R(h) 



We have finally obtained a differential equation linking strain (contained 
in h = (utk — 5)/2) and the crack tip position xq. The equation is nonlinear, 
and we cannot give an analytical solution in a general case, but many general 



3 Same as 1 
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features can be worked out. First of all it can be checked that the right hand 
side of Eq. ( 121 j) vanishes linearly at 8 = 8q. This indicates that the crack 
tip advance X (8) has a logarithmic divergence at 8q. Second, we note that 
Eq. (l2Tj) will determine the crack tip position as a function of 8 up to an 
additive constant, i.e, the solution will be of the form 

x = F(8) + A, (22) 

with F(S) a well defined function diverging logarithmically at 8q , and A an 
arbitrary constant. The value of A can be set by using the condition that X 
becomes different from zero exactly at 8q for a virgin sample. This defines 
the solution unambiguously. 

Expanding R(h) (see equation (FIT]) ) in powers of h, we obtain to linear 
order 



which becomes exact as (P — > 0). 

Using this approximation for R(h) in (12 ip leads to the following linear 
first order differential equation 



^ = -(^ + VELU + "» (24) 

In this approximate case, the solution is fully described by a logarithmic 
dependence of xq as a function of 8, rising linearly from zero at 8q, and 
diverging at 8q. 

To check this behavior, a comparison of the analytical expressions ob- 
tained with results of numerical simulations is presented in Fig. |HJ The 
numerical simulations show an important dependence on t he parameter A 



which can be considered in terms of a lattice trapping effect (IThomson et al 



1971c iPaskin et al.1 . Il98ll ; iKessler and Levind . Il999al ; iBernstein and Hessl . 120031 ) 



The continuous limit is obtained by letting A — > 0. We see that when this 
is taken into account the agreement between the analytical and numerical 
results is very good. 

In cases in which we cannot use the linear approximation (1231) for R(h), 
Eq. (]2ip has to be treated numerically. For instance, the slope of the right 
hand side of Eq. (]2ip at zero crossing determines the pre-factor of the log- 
arithmic divergence at 8q . The comparison of this pre-factor with results 
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Figure 6: (a) Analytical results for the crack tip advance X according to the linear 
approximation given in Eq. (|24[) and results of numerical simulations at different values 
of A. (b) The same results shifted according to the values of 8q (v,bk/u n i = 4, P = 0.5). 
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obtained from full numerical simulations of the model are contained in Fig. 
([7j). We observe an almost perfect agreement when using Eq. (12 ip . and 
deviations if using the linear approximation Eq. ( 124|) . 
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Figure 7: Results for the pre-factor a of the logarithmic divergence of X(5) at 5jJ . Circles 
are obtained fitting exponential functions to our numerical results with Ubk/u n i = 4, and 
P = 0.5. Continuous line is obtained from the full Eq. (f2Tj) . and dotted line is the result 
of the linear approximation (Eq. [24]) . 

3.2. Cyclic crack growth and Paris curves 

Crack tip advance upon monotonous increase of the applied load, anal- 
ysed in the previous section, is the starting point for the study of the more 
interesting regime of cyclic loading conditions. Particularly we consider the 
case in which, after the initial application of some S max , the stretching is 
reduced to some 5 min , and then it is successively cycled between these two 
values. We expect a finite crack propagation in each cycle of this process. 
For this to happen, a crucial condition must be fulfilled: in the unloading 
part of the cycle, some of the plastic deformation must be reverted. Oth- 
erwise, if downloading is totally elastic, the next loading half cycle will be 
elastic too, and crack tip advance will not occur. Reverted plasticity allows 
a sustained advance in successive cycles. In fact, imagine the hypothetical 
situation in which all plastic deformation induced during the first loading is 
reverted upon unloading, namely the Iq of all springs reset to zero value. This 
would mean that after a complete cycle, the system is back in the original 
"virgin" state, with the only difference that crack tip has moved forward, 
and the situation will be repeated cycle after cycle. 
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The numerical results within our model fit in between these two limits. 
They are presented in Fig. [SJ Starting from a virgin sample, we first reach 
some value 5 ma x in between 6q and Sq. This generates a first crack elongation 
and a plasticity wake that was discussed in detail in the previous section. 
Then 5 is reduced to 5 m i n . If during this reduction of 5 some part of the 
plasticity in the system is reverted, the system can elongate the crack further 
in a second increase of 5. Whether this plasticity reversion occurs or not, 
depends a great deal on the value of b~ m in- If this value is too high it does not 
occur. In the case of the present model, we have seen that we have to take 
b m i n negative in order to have reverted plasticity. For this reason in Fig. [8] 
we use 5 min = -5^. 

A brief digression is convenient at this point. The use of 5 values of both 
signs is perfectly allowed in our mode III configuration. However, having 
in mind an experimental situation in mode I configuration, where 5 should 
be strictly possitive, we may worry about the neccessity to include values of 
S of alternating signs in the present case in order to observe fatigue crack 
propagation. In this respect, we mention that the need of a change of sign 
of S is an artifact of the quasi one dimensional system. In two dimensional 
(i.e., many chains) mode III simulations we have observed reverted plasticity 
and cyclic fatigue crack advance for strictly positive values of 5 m i n . 

In Fig. [8] we see that part of the plasticity is reverted during S reduction. 
We stress, however, that for the current parameters there is no crack advance 
during the stress reduction half period. Then upon a new increase of 5, crack 
elongation resumes at some 5 < 5 max . The form of the xo(S) crack tip advance 
curve in successive cycles can be described by assuming that previous cycles 
only influence the advance in the present cycle by shifting the value of 5 at 
which elongation resumes. This is understood in terms of the analysis of 
the previous section: Eq. ( 12Tj) is still valid to analyze the crack tip advance 
in the presence of plasticity behind the crack tip. The only difference is 
that now the determination of the constant A in Eq. (122|) cannot be done a 
priori, and the plasticity in all the region behind the crack tip has to be taken 
into account. However, once the crack tip starts to elongate, its evolution is 
dictated by Eq. (I2T|) . We will not attempt here to explicitly calculate the 
value of 5 at which elongation starts in each cycle, in terms of the plasticity 
distribution behind the crack tip. Instead, we restrict to a rather qualitative 
description. 

After a few cycles, the elongation of the crack converges to a fixed amount 
per cycle. The asymptotic value of advance per cycle is seen to be lower than 
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the advance in the first cycle, when there is no plasticity in the system. 
This is consistent with the fact that plasticity is only partially reverted upon 
unloading. Serrations in the profile of plastic deformation are clearly visible 
in lo(x) (Fig. [8]). They manifest also in the chain profile itself, as the last 
panel of the figure allows to observe. This is reminiscent of the same texture 
observ ed in crack surf aces generated by cyclic fatigue (see for instance ISuresh 
fll998h . Fig. 10.4, and iHertzberel ( Il996h . Fig. 13.11). 

By systematically running the model at different values of 5 m i n and 5 max , 
Paris-like curves can be constructed, displaying the stationary advance-per- 
cycle as a function of load amplitude. In Fig. fl9]) we present these results 
as a function of S max , for four different minimum load values £ m j n , in the 
usual Paris form, and also in a more appropriate display to our analysis. We 
observe that the full curves can be reasonably described by a logarithmic 
dependence with Sq — S max , with only a rigid vertical shift to account for 
the different values of S m i n , namely, the advance per cycle dx/dn has the 
qualitative form 



dx 

— ~ a \og{S, 



ii 



+ m 



(25) 



This expression cannot be put in the form of a dependence on a single 
combination between S max and S m i n , as traditional interpretations of the 
Paris law require. It is however compatible with two-parameter interpre - 
tations of the Par i s law, as suggested for instance by Vasudeyen et al.l (| 19941): 
Sadananda et al.l (ll999l ): ISadananda and Vasudevanl (120041 ) : IVasudevan and Sadananda 

(Hqo3). 

Note that there is a well defined value of 5 max (depending on S min ) below 
which there is no systematic crack advance. This indicates the existence of an 
endurance limit for our model below which fatigue crack propagation does not 
occur. The rate of crack advance increases with 5 max , and becomes very large 
approaching Sq . In the middle, there is a transition between the two limiting 
cases. However, this intermediate regime (which is the typical Paris regime) 
is very narro w compared to most experimental situations (see for instance 
Broekl (119821 ). Fig 10.1). We associate this fact to the existence of a single 
chain in our model. If the system was "two dimensional" (i.e., with many 
chains), plastic deformation would appear in all the chains, and the load 
region in which fatigue crack advance occurs would be largely increased. We 
have obtained this enhancement effect in two dimensional simulations, and 
plan to report on it elsewhere. For our present simulations, it is remarkable 
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that such a simple model displays the same qualitative behavior observed in 
a large class of real materials. 



3.3. Overload retardation 

A characteristic effect observ ed in cyc l ic cra ck growth is overload retarda- 



tion (Sureshl ( Il998h . pp520-526, iBolotin! dl999|)pp.257-263, ISadananda et al 



( 119991 )). Consider a crack subjected to cyclic loading, advancing some finite 
distance per cycle. If the maximum load in one particular cycle is increased, 
the crack will advance a larger distance in that cycle. But also a relatively 
bigger plastic deformation will be produced, that will shield the crack tip 
and reduce the advance in the following cycles so that eventually, when the 
system reaches steady state advance again, the crack tip position may be 
lagged with respect to the case in which overload was not present. In some 
cases, the crack may even completely arrest after the overload. This some- 
what counter intuitive phenomenon is very well known experimentally, and 
is a stringent constraint for any model that is supposed to describe fatigue 
crack growth. Models that have been proposed to describe this behavior 
may be c l assified into two main catego ries. In crack tip plasticity models 
( jWheelerl . Il972t IWillenborg et al.l . Il97ll ) it is assumed that crack growth re- 
tardation occurs due to the large plastic zone developed during overloading. 
The residual compressive stresses formed in this zone will reduce the mag- 
nitude of the tensile stresses during the next fatigue cycle and tend to delay 
crack growth. In crack closure models (jElberl . Il97ll ) it is argued that as a 
result of the tensile plastic deformation left in the wake of a fatigue crack, a 
partial closure of the crack faces occurs during part of a fatigue load cycle. 

We have been able to reproduce the overload retardation effect using our 
model (Fig. (fTUj) ). The system is initially cycled between S max = 1.0545^ 
and 5 m i n = —0.65q. After about n « 10 cycles the growth rate stabilizes at 
dx/dn = 0.2 which corresponds to four lattice units per cycle, taking into 
account that A = 1/20 in this simulation. In cycle n = 11 the overloading 
is applied, through an increase in S max of about 1%, and the crack is seen to 
advance more than 3 times the previous value. In the following cycles the 
crack advance is seen to be reduced drastically, and the whole effect once the 
original rate is recovered results in a net retardation of 2 lattice units. The 
reason for this reduction in crack tip advance is the larger plastic deformation 
induced during the overload cycle. Fig. f flUjl b shows this enhanced plastic 
deformation in lo(x) at the overload cycle. The additional plastic deformation 
has a stabilizing effect on the crack tip, so the application of 5 max after the 
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Figure 8: (a) Temporal evolution of crack tip position upon cyclic change of load, starting 
with a virgin sample. We see that after a few cycles, the advance per cycle stabilizes to 
a finite quantity, (b) Snapshots of the system during the process. Panels are snapshots 
at the first d m ax, the first 5 m in, the second S m ax, and at 5 max after many cycles. The 
typical serrations of cyclic fatigue advance are observed in the plastic deformation profile 
behind the crack tip and in the geometrical profile of the chain itself, as panel (c) shows. 
Parameters used are P = 0.5, A = 1/20, Ubk/u„i = 4, 5 V 



= — Sq and 8„ 
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Figure 9: (a) Stationary crack tip advance per cycle dx/dn as a function of 5 max , for 
different values of S m i n , as indicated. The curves show an overall form qualitatively com- 
patible with a Paris law. In (b), we see that results for different values of 5 m i n can be 
absorbed by a vertical shift of the different curves (the vertical shift applied is indicated). 
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overload is much less efficient in generating crack tip advance. Eventually, 
when the crack tip has moved far away of the enhanced plasticity zone, the 
crack tip elongation rate returns to its equilibrium value. 




Figure 10: Overload effect: a) crack tip position versus number of cycles. The crack is 
loaded cyclically between 5 max — 1.054<5<t. and 8 m i n = —0.66q. At n = 11 an overload 
^overload = 1.063(Jg is applied. The effect after many cycles amounts to a net retardation 
of two lattice units, in the present case, b) Profile of the system u(x) and the springs rest 
length Iq(x) in the final configuration (simulations with A = 1/20). Compare panels a) 



and b) with Figures 10.10, and 2.29 in lBroekl (|1982l ). 



Whether the combined effect of the larger advance in the overload cycle 
and reduced advance in ulterior cycles gives an overall retardation or not, 
depends in a delicate way on the parameters and the values of 5 that have 
been used. However, one extreme case in which we can be more quantitative 
is the following. We investigated what the conditions for the largest possible 
overload to arrest the crack advance are. To answer this question, we first 
refer to Fig. [TTJ There, the continuous line (picked up from data as that in 
Fig. [9]) separates the region of cyclic crack propagation (above the curve) 
from that of no propagation (below). Now, on each situation of a propagating 
crack, we apply the largest overload the system is able to sustain, namely, one 
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in which 5 is increased to almost 5q @. This produces a very large elongation 
during the overload cycle. After that, 5 is reduced to the corresponding 
S m i n , and then increased. The value of 5 at which crack tip starts elongation 
again (indicated by the dotted line in the figure) is the minimum value of 
Smax required for crack elongation after the overload. In other words, in the 
region between continuous and dotted lines in Fig. [UJ a crack is arrested by 
the overload in the system. This quantifies appropriately the conditions for 
crack arrest in the system, showing clearly that conditions can be found for 
the process to occur. Typically, above the curve of crack arrest, there is a 
region of overload retardation (shown shadowed in Fig. [TTj) that we have not 
determined in detail for all the parameters. For even larger values of S max , 
a net retardation of the crack does not occur, instead the crack gets a net 
advance due to the overload. 




-1.6 -1.4 -1.2 -1 -0.8 -0.6 -0.4 -0.2 



Figure 11: The plane 5 max vs S m i n , indicating the fatigue limit (continuous line) separat- 
ing regions in which cyclic crack propagation occurs, or not, for the current parameters 
(ubk/uni = 4, P = 0.5). Within the region of fatigue propagation, the dotted line limits the 
region in which this propagation can be arrested by a single overload in the system. Above 
this line, a region of overload retardation (dashed, not accurately determined) exists. 



4 Note that as this overload generates a very large elongation, the result is independent 
of previous values of S m i n and 8 max . This means that a single simmulation of the overload 
sufhcces for all 8 m i n and 8 max . 
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4. Conclusions 

In this paper we have introduced a simple quasi-one-dimensional model 
of a crack advancing in between two elasto-plastic chains. The chains are 
joined by breakable springs and attached to two lateral strain gauges by which 
external loading in the form of prescribed strains are applied to the system. 
The use of elasto-plastic elements forming the chains gives the possibility of a 
cyclic fatigue crack advance, in which a finite advance of the crack is obtained 
at each cycle of the external strain between a maximum and a minimum 
value. For the case in which we start with a virgin sample (i.e., without 
plasticity) we have derived analytical expressions to solve the position of the 
crack tip in terms of the applied strain, and have successfully compared these 
results with those of numerical simulations. In the case of cyclic advance, 
we have numerically obtained curves for the amount of elongation in terms 
of the maximum and minimum strain. We have observed the existence of 
an endurance limit below which fatigue propagation does not occur, and a 
maximum strain above which the abrupt rupture of the sample occurs. 

The mechanism of fatigue crack advance in the present model is crack tip 
shielding by a plasticity wake behind the crack tip. It can be described in 
the following way. During load increase, a stable crack tip advance occurs 
due to the generation of a plasticity wake behind the crack tip that avoids 
the immediate unstable propagation of the crack. Upon load reduction, part 
of the plasticity in the vicinity of the crack tip is reverted. In the second 
period of load increase, the zone of reverted plasticity allows for a new finite 
elongation of the crack. As the reverted region is typically smaller than the 
plastic region induced during crack advance, the elongation slightly decreases 
in consecutive cycles until it reaches a stationary value. The plasticity profile 
left behind by the advance of the crack has typical serrations that reflect 
in corresponding undulations of the chain profile, reminiscent of the same 
phenomenon during cyclic fatigue propagation in real materials. 

Crack tip advance per cycle was seen to be a function of both the maxi- 
mum and minimum value of applied strain 5 max and 5 m i n . This means that 
it cannot be simply written in term of a single AK, as in the simplest form 
of a Paris law. It could be interpreted that an additional dependence on the 
stress ratio R = K min /K max exists. If a single parameter interpretation of 
the fatigue advance can be given in our case, this single parameter should be 
considered to be the amount of plasticity reverted during the unloading part 
of the cycle. However, this amount depends on both S max and 5 m i n giving a 



25 



two-parameter Paris law. 

In addition, we have studied and reproduced the effect of overloading, 
showing that a cycle with an excess of applied strain can lead to a net retar- 
dation in the elongation of the crack, and even to a complete arrest of it. We 
have numerically determined the conditions in the plane S max ,5 m i n for this 
arrest to occur. The origin of crack retardation or crack arrest in our model 
is the shielding of the crack tip by the enhanced plastic deformation during 
the overload cycle. 

The minimal model we have presented allows a detailed analysis of the 
fundamental processes that lead to the possibility of fatigue crack propaga- 
tion. In a forthcoming publication we plan to discuss the similar properties of 
a more realistic system consisting of a two-dimensional mesh of elasto-plastic 
springs and masses loaded in a mode I configuration, where qualitatively the 
same kind of fatigue propagation is observed. 
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